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Abstract: We study Maxwell-Chern-Simons theory in 2 noncommutative spatial dimensions 
and 1 temporal dimension. We consider a finite matrix model obtained by adding a linear 
boundary field which takes into account boundary fluctuations. The pure Chern-Simons has 
already been shown to be equivalent to the Laughlin description of the quantum Hall effect 
1^, 0]. With the addition of the Maxwell term, we find that there exists a rich spectrum 
of excitations including solitons |^] with nontrivial "magnetic flux" and quasi-holes with 
nontrivial "charges", which we describe in this article. The magnetic flux corresponds to 
vorticity in the fluid fluctuations while the charges correspond to sources of fluid fluctuations. 
We find that the quasi-hole solutions exhibit a gap in the spectrum of allowed charge. 
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1. Introduction 

In recent years noncommutative geometry has been an effervescent field of research particu- 
larly in its relation to solitons in effective descriptions of string theory and D-branes 
However its most surprising application comes in a description of strongly correlated quantum 
magneto- hydrodynamics and various other quantum dynamical fluids |^. The most intrigu- 
ing application in this context is to the quantum Hall effect. SusskindQ proposed that non- 
commutative Chern-Simons theory in 2+1 dimensions would be an appropriate description of 
the quantum Hall effect. The quantum Hall effect concerns the strongly correlated quantum 
dynamics of a two dimensional electron gas in a strong transverse magnetic field. The non- 
commutative space exists in the internal two dimensional space of the Lagrange coordinate 
description ||5| of the electron fluid. 

The continuum, classical description of the small fluctuations of a two dimensional fluid 
is easily seen to be a gauge theory of the group of area preserving diffeomorphisms. The gauge 
fields (spatial components) correspond to fluctuations of the fluid with respect to the ground 
state of a quiescent, undisturbed fluid. The gauge freedom of area preserving diffeomorphisms, 
simply corresponds to a relabeling of the elements of the fluid which are native to the La- 
grange description of fluid dynamics, an evident invariance of the theory. The corresponding 
conservation law is equivalent to the Gauss law. 

In the presence of a strong transverse magnetic field and in the low energy approximation 
the classical term with the lowest number of derivatives is exactly the Chern-Simons term. In 
this theory, the Gauss law in fact imposes the vortex free condition on the fluid. The vortices 
are frozen out of the fluid and act as sources, just as ordinary charges in electrodynamics 
act as sources outside of the electric and magnetic fields. Imposing the Gauss law via the 
introduction of a temporal gauge fields and an enhanced gauge invariance (now including 
time dependent gauge transformations) results in a fully non-abelian Chern-Simons gauge 
theory of the group of area-preserving diffeomorphisms. 

Susskind's Q key observation was that this non-abelian gauge theory appears to be a 
truncation to first order of the simplest noncommutative Chern-Simons gauge theory defined 
on two noncommutative spatial and one normal temporal dimension. Thus Susskind pro- 
posed that the true theory of the quantum Hall effect corresponds to the full noncommutative 
Chern-Simons gauge theory. One motivation given for this hypothesis was to reintroduce the 
discreteness that exists at the particulate level of the two dimensional electron gas, a discrete- 
ness which the continuum approximation erases. It remained to be seen if the phenomenology 
of the quantum Hall effect could be reproduced with this hypothesis. It was shown in Susskind 
but also in more detail in ||6|, ^] that indeed the noncommutative gauge theoretical descrip- 
tion of the quantum Hall effect when restricted to a finite droplet of the fluid through the 
introduction of a boundary and boundary degrees of freedom, was in one to one correspon- 
dence to the description afforded by the Laughlin wave functions Q. However, the probability 
densities calculated in the noncommutative Chern-Simons model was only equal to that of 
the Laughlin wave functions in large distance limit ||9|. This theory however describes the 
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quantum Hall state via the projection to the lowest Landau level, it cannot hope to describe 
any transition between levels or the transition to the final state called the Hall insulator |1C] 
for very strong external magnetic field. 

In the absence of the transverse magnetic field, the lowest order term in the effective 
Lagrangian corresponds to the Maxwell term for the gauge field of area-preserving gauge 
transformations. The coefficient of the "electric" part does not have to be correlated with the 
coefficient of the "magnetic" part combining to give a relativistically invariant action, however 
this can be arranged by re-scaling the gauge field or time variable appropriately. Again the 
Gauss law constraint can be obtained by the incorporation of a temporal component to the 
gauge field and writing a gauge invariant expression for the field strength. The Maxwell 
term is the next order term that can be added to the pure Chern-Simons gauge theory. It 
renders the theory more interesting, the Gauss law constraint does not expel vortices from 
the theory but imposes a more dynamical constraint. We have already studied this theory 
in previous articles, where we found plane wave solutions for the unbounded theory |11| and 
soliton solutions for the theory of a finite droplet |12|. In this article we further examine the 
theory on a finite droplet and show the existence of quasi-hole states and rotational excitations. 
With this rich spectrum of excitations we expect that the theory should be able to describe 
transitions as a function of the parameters. Other authors|jl3| have looked for the quantum 
solution of the noncommutative Maxwell-Chern-Simons theory and found the correspondance 
to be to more than one Landau level. However, in order to find their solution they had to 
assume that a certain deformation energy of the fiuid (defined in the next section) was either 
zero or infinity. In this article, we find classical quasi-hole solutions for the noncommutative 
Maxwell-Chern-Simons theory for arbitrary values of the deformation energy. 



2. The model, equations of motion and the Hamiltonian 

Susskind'sQ idea was to describe a two (spatial) dimensional fiuid by a gauge field Aj so that 

x^ = y' + ee'^Aj (2.1) 

where Xj, i = 1,2 are the Eulerian coordinates of the fiuid, yi, i = 1,2 are the Lagrangian 
(comoving) coordinates of the fiuid and 9 = l/{2'Kpo). Then the Lagrangian of a charged fluid 
in an external transverse magnetic field corresponds to the 2-1-1 dimensional Maxwell-Chern- 
Simons theory for small value of the gauge field. The continuum approximation removes the 
discreteness that is manifest in the physical fiuid. Susskind proposed to bring this discreteness 
back by suggesting that the noncommutative version of this theory should describe the full 
theory. We studied this theory with the additional modification of boundary and boundary 
degrees of freedom[^]. Here we study the same action, however we add a factor k to the 
analog of the magnetic field squared term which corresponds to the potential energy density 
of spatial deformations of the fiuid. The first term, the analog of the electric field squared, 
corresponds to the kinetic term, the Chern-Simons term represents the interaction of the 
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charged fluid with the external magnetic and the last term represents the boundary degrees 
of freedom: 

^ = ^/ dt{Tr{{-2[D^,D][D^,D^]-K[D,D^][D,D^])) 
+2X{-[D,D^] + 1)Dq] -2^^Dq^). (2.2) 

Here Dq is the time covariante derivative, D and D'^ are the holomorphic and anti-holomorphic 
combinations of the spatial covariant derivatives respectively and ^' is a boundary field. The 
boundary field was first added by Polychronakos ||6| which allowed him to find solutions of 
non-commutative Chern-Simons theory in terms of finite matrices. These correspond to finite 
droplets of the quantum Hall fluid. Specifically 

= ^fe{-i^^ + A,,) (2.3) 

°=^-"' = ^ (2.4) 
and the parameters A and are related to the noncommutativity parameter 0, the electron 
charge e, the external magnetic field the density pQ and the electron mass m by 

A=£^, = (2.5) 

m m 

We rescale Dq and ^ and the parameters in the following way, in order to obtain exactly the 
action studied in IIT2|: 



K mJK rriK \/k 



Defining 

we obtain the action 



(2.7) 



S = -Ej dt{Tr{{-2[D^,D][Do,D^]-[D,D^][D,D^])) 
+2X{-[D,D^] + 1)Dq] -2^^Dq^). (2.8) 
By varying with respect to Dq, and we get, respectively, the boundary equation 

iir = Aq^, (2.9) 

the Gauss law 

[D, [Do, Z)t]] + [D\ [Do, D]] + X{[D, D^] - 1) + *M/t = 0, (2.10) 
and the Ampere law 

[Do,[Do,D]] + [D,[D,D^]] = X[Do,D]. (2.11) 
The Hamiltonian is, as in ||T^, given by 

H = E Tr{-2[Do,D^[Do,D] + [D, D^][D, D^]). (2.12) 
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3. Rotational excitations 



Our first solution corresponds to rotational excitations on top of any given solution. Our 
procedure can be applied to the soliton solutions found for example in |12| and to the solutions 
that we find in this article. We put 

R 



D = D' + 



Awt 



26 



Da 



D'o, ^ 



(3.1) 



where the primed variables correspond to any known solution to the equations of motion 
with R (proportional to the identity) and w (real) simply constant. The Gauss law involves 
commutators of the D or D^ with each other or with their commutator with Dq. The direct 
addition of complex constants as in the equation (3J.) or those that result from the commu- 
tators involving Dq, simply vanish open taking the further commutators with D or D\ hence 
the Gauss law is satisfied. The equation ( ^^ ) is also obviously satisfied. Replacing (3^) into 
the Ampere law yields a solution ii w = = ^ the familiar cyclotron frequency. We 
calculate the corresponding change of energy: 



R* 



AH = -2ETr [Do, -=e~™*] [Do, I?'] + [Do,D'^[Do, 



R 



29 



+ [^0, 



R* 



-iwt 



29 



R 



Awt 



which yields 

AH 



2ETr ( -xJ^e-''^^[Do, D'] + A 



R 



e-[Do,D'^] 



Awt r 



A- 



,R*R 



26 
(3.2) 



(3.3) 



Since the commutator [Dq,D\ is off-diagonal in our solution and the solution in ||12| the trace 
vanishes giving 

^=^'1^1^'^. (3.4) 



AH 



29 



Expressing this in terms of physical constants gives 

e'^B^[R['^N 



AH 



2m 



(3.5) 



which is exactly the energy of N rotating electrons with a amplitude R at the frequency w in 
a magnetic field. 

This solution adds a term to the covariant derivative 



D[ + Li cos{wt + (p),D2 = D2 + sm{wt + (p) 



Di = D'l + ^-j^ cos{wt + ip),D2 = D2 + sm(wt + 93) (3.6) 
where f is the phase of R. Using the Susskind correspondence between the fluid coordinates 



and the gauge field 
original solution as 



we see that the solution (3^1) corresponds to rotating about the 



yi=y[-[R[ sm{wt + (p), y2 = y'^ + [R[ cos{wt + if) (3.7) 

where y[ and 2/2 are for operator coordinates of original solution. These rotational excitations 
clearly also exist in the analogous theory in the infinite plane treated in |11|. 



4. Quasi-hole solutions 



Quasi-hole solutions appears as a modification of the solutions found in |12|. We will take 
an ansatz similar to |12|, hence D is represented by an x matrix which satisfies certain 
boundary conditions, however we will take a periodic D as in [||]. The D is an operator similar 
in structure to the annihilation operator of an ordinary Heisenberg algebra, hence it generally 
relates a state |n + l) to a state |n). For the finite matrix representations used here, periodicity 
means that the final state |0) is related to the state — 1) since n ranges over the N values 
0,1,--- A^ - 1. 

D=2Z ^G{n) +qe Ve | n) (n + 1 | +^e^^ | iV - 1) (0 | (4.1) 

n=0 



If g = 0, this ansatz is equivalent to the one in |12|. G(n) is to be determined from the 
equations of motion, we solve these eventually, perturbatively and numerically. 

Solutions with this ansatz correspond to quasi-hole solutions because they bound the 
lowest eigenvalue of the radius away from zero in a g dependent manner. From (|2.1|), the 
fluctuation of the radius is proportional to the square of the gauge field 

Ai,Al.{^^,.0.)\[^^,,^;)\ ,4.2, 

Noncommutative geometry in a finite space is defined by the commutator of the coordinates 
[xi, X2] = i0{l — N \ N — 1){N — 1 I) which imply the commutator for the derivate [di, 82] = 
|(-1 + N \N -1){N Thus we can put 

ft = -i=(. + .t),4 ,4.3) 



where 



Then we obtain 



N-2 



d = ^ Vn + 1 I n) (n + 1 I . (4.4) 



n=0 



oc{D - d,D^ - (4.5) 
where is the anti-commutator. Then the radius becomes 



oc (g{0) + 2g - 1 - 2y/G{0) + gcos | 0)(0 

J 



n=l 



+ (^G{N -2)+2q-N + l- 2^G{N -2)+q^W^cos { ^^^ ^^^^ ) ) I ^ " 1>(^ " 1 

(4.6) 
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a diagonal expression in the states where A{n) = A{n) + A{n — 1). We see that for large q we 
have correspondingly large eigenvalues for R^. The smallest eigenvalue is not directly equal 
to q, hence our solution corresponds to the fluid pushed away from the origin. 

Returning to the solution of the equations of motion, we will consider the gauge = 
but we let D depend on time. This choice is different from that taken in |12|, however in 
that case, our choice is simply gauge equivalent. With the periodic ansatz of equation (4.1), 
this is not the case. Indeed, our choice give us an additional degree of freedom which we can 
identify as p. Thus Dq = —iVOdo, and we can calculate the different terms in the equations 
of motion: 

N—2 

[Do,!?] = ^ ^i;(n)VG(n)+ge ^ | n) (n + 1 | +p^e |7V-1)(0| (4.7) 

n=0 
N-2 

[Do, [Do, D]] = ^ (u;(n)) VG'(n) + qe Ve \ n){n + 1 [ ^ | - 1)(0 | (4.8) 

n.=0 

The commutator [D, D^ is given by replacing from equation (|4.l|) 

[D, D^] = 22Y1 V^^(") + iVGim) + qe ^ [| n)(n + 1 |, | m + l)(m |] 

n=0 m=0 

+g[|iV-l)(0|,|0)(iV-l|] 

N-2 N-1 

= ^(G(n) +g) I n)(n | - ^(G(n- 1) +g) | n){n | +q{[ N-1){N -1 | - | 0)(0 |) 

n=0 n=l 
N-l 

= J2G{n-l)[ n){n [ (4.9) 

n=0 



with the notation A{n) = A{n + 1) - A{n) and we define G{-1) = and G{N -1) = 0. The 
commutator [D^ , [Do,D]] can be computed in the same way as [D,D^. Defining u;(-l) = p 
and w{N — 1) = p then 



N-l 



[D^,[Do,D]] = - ^ {G{n - 1) + q)w{n - 1) [n){n[ . (4.10) 



n=0 



As this result is hermitian it is also equal to [D, [Do,D^]]. The final term appearing in the 
equations of motion is calculated as 



N-2 N-l . , ^ ^, 



[D, [D, D^] = ^Y1 y^G{n) + qe ^ G{m - l)[[ n){n + I [,[ m){m [] 

n=0 m=0 
^^^^ N-l 

+^e^ ^ G{m -l)[[N - 1)(0 |, | m){m |] 
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N~2 



VG{n) + qe 



iw ( n) - 

Ve G{n 



1) I n)(n + 1 I ^ 



n=0 



G(-l) -G(7V-2) I iV 



1)(0| 
(4.11) 



with the notation G{n -\) = V'^G{n) = G{n + 1) - 2G{n) + G(n - 1) where V'^G{n) is 
discrete Laplacian (which is facihtated with the further notational convenience G{N) = G(0) 
and (A^ |= (0 |). Then 

[D, [D,D^]] = Y VG{n) + qe ^ V^G{n) \ n){n + l\ . (4.12) 
?i=0 

The solution for ^' is simply a general static vector since = in ( ^.91 ): 

Af-l 

* = ^ A„ I n). (4.13) 

n=0 

We see, as in |12|, that only contributes off diagonal terms in the Gauss law ( |2.1ClD . 
To eliminate such terms we must take ^ = Xm \ M). Contrary to |12|, different choices of 
M are all gauge equivalent (i.e. a permutation) since all choices of M are equivalent in our 
periodic ansatz (up to the name of the variable). Taking the trace of the Gauss law then 
yields Am = VNX. We choose without loss of generality that M = N — 1. Then the Gauss 
law yields (for n = 0, • • • — 1) 



-2(G(n - 1) + q)w{n - 1) + A [G{n - 1) - ij + XN6N-i,n = 0. 

We will solve this equation by induction as in jl^. We will show that for n = [0, A'' 

n + l-^P- G{n) 



w{n) 



-X- 



(4.14) 



(4.15) 



2{G{n)+q) ■ 

This formula is true for n = — 2, as is verified by considering the Gauss law ( 4.14| ) for 
n = A^ — 1. Then assuming the form (fl.lSj) for a general value of n we can prove that it 
is valid for n — > n — 1. Thus by the principle of induction the formula is valid for all n. 
However we have to check/impose that the n = equation of ( [4.14| ) is respected. This is 
indeed the case if we take ( [4.15| ). Taking the Ampere law ( p.lip and removing an overall 



factor Y^G(n) + ge we obtain, for n = [0, A^ — 1]. 

{w{n)f +V'^G{n) = Xw{n). 



(4.16) 



Here we have A^ equations, in A^ — 1 values of G{n) and also in the two variables q and p. Hence 
we have A^ equations and A^ + 1 parameters, thus we have one free parameter. Generically, 
there will be a family of solutions. The Hamiltonian, from equation 2.12 , for this ansatz, using 
equations 4.7 and 4.9 is 



N-l 



H = E{Y 2{w{n)f{G{n) + q) + (G{n - 1)) 



(4.17) 



n=0 
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4.1 Perturbative analysis of the quasi-hole solution 
4.1.1 Large quasi-hole solution 

The equation ( 4.16| ) is non-linear and thus difficult to solve. So we will look the solution for 
g 3> 1. We define 

p = po + ^, G{n) = Go{n) + (4.18) 

q q 

Then if we expand w to first order in 1/q, we obtain: 

win) ~ po - T^(A(n + 1 - Go(n)) - 2pi + 2/5oGo(n)) (4.19) 
/q 

{w{n)f ~ p2 _ ^(^(^ ^ 1 _ _ ^ 2poGo{n)) (4.20) 



Working to zero order in 1/q, equation 4.16 gives (for n = [0, — 2]) 

/Oo + V'Go(n) = Apo- (4.21) 
The discrete Laplacian solved as in the continuous case by 

Go(n) = aon^ + p^n + 6o, ao = (4.22) 

We have to impose the two boundary conditions Go(— 1) = Go{N — 1) = 0. These give 

ao-Po + 6o = 0, ao{N - if + Po{N - 1) + 5o = (4.23) 

which are solved by 

/3o = ^(iV - 2)(p2 _ Apo), 5o = ^(7V - l)(pg - Apo). (4.24) 
We still have one final condition left (from the n = — 1 of equation 4.16| ) which gives 

pI + ao{N - 2)2 + (3o{N - 2) + 26o = Apo (4.25) 

with solutions 

po = or po = A. (4.26) 

Either of these solutions give us Go{n) = 0. 

With the order zero solutions we can continue to solve the equation ( [4.16| ) in first order 
in 1/g 

-po(A(n + 1) - 2pi) + V^Giin) = ^(A(n + 1) - 2pi). (4.27) 

This is solved by 

Gi(n) = am^ + /3in2 + dm + 71, ai = ^(2po - A), /3i = ^(2po - A)(A - 2pi). (4.28) 
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Again we have our boundary conditions 1) = Gi{N — 1) = which imply 



Si = ^iN''X-3X- mpi + 12pi) (A - 2/5o) , 71 = - ^ (A' - 1) (AiV - 6/)i + A) (2po - A) . (4.29) 



Finally the condition for n = N — \ of ( [4.16 ) gives 



Thus 



Pi = ^(iV-l)A. (4.30) 



'^1 = -^(A^'-9A^ + 12)(A-2po) (4.31) 
71 = ^(iV-l)(iV-5)(2po-A) (4.32) 
/3i = ^(3-iV)(2po- A) (4.33) 



From the Hamiltonian ( 4.17| ), for the po = solution, it can be shown with a little calculation 



that a non-zero contribution arises only at order 1/q. The variables {p,G{n),w{n)) must 
be expanded to order 1/g^ to consistently extract this contribution, because the Hamiltonian 
contains an explicit factor of q. Indeed, as we will see below, for the po = solution, the second 
order expansion of these variables does not in fact give a non-zero contribution, however they 
do contribute to the energy for the po = A solution. Then we find the solution to second order: 

G2(n) = a2n^ + p2n^ + S2n^ + 72?^^ + 6?^ + £2, (4.34) 
a ^ (3 ^ -a + 25 
= 20' = I^' = -^T- 

6 = + ^(3a - /3) + —(-5a + 4/3 - 25) (4.35) 

N 1 
+ —(3a - 5/3 + 65 - 67) + -(/3 - 35 + 67) 
12 6 

-/? + 67 , g - 55 + 157 
72 = 62 = 6 , (4.36) 

a = /? = ^(2/3i - 1), 5 = ^(2pi + A5i - A) 

7 = Api - + Ap2 - y - 2p2Po + ^ (4.37) 

We have implicitly assumed that q ^ \, q ^ pi and q » Gi{n) for every i (or more precisely 
that these variables are of order of q^). However, we find that the Gi{n) are polynomial in 
n where n = 0, • • • , A^. Thus if N becomes large, Gi{n) would also become large. Therefore 
the condition g 1 is not adequate for the perturbative expansion to converge. We see that 
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Gi{n) and G2{n) are respectively third and fifth degree polynomials. Thus our perturbation 
series would not be valid if A^'^ 3> g or A^^ ^ q^. If we replace 



G{n) = Y^^. (4.39) 



i=0 



in equation ( 4.16| ) and expand to order l/g*, then we obtain a recurrence equation relating 



the discrete Laplacian of Gi{n) to Gj{n) with j < i. Assuming that the solution for the Gj{n) 
with J < z is a polynomial in n, the solution for Gi{n) is a polynomial for degree two higher. 
Then knowing that Gi{n) is a polynomial of order 3, by induction we see that Gi{n) is a 
polynomial of the degree 2i + \. Thus we see that the true perturbative parameter is actually 
is -n? /q. Since n = 0, • • • , A^, we take the strongest condition, N"^ /q <^ 1. This condition might 
actually not be necessary, since we in fact impose the boundary condition that G{N — 1) = 0, 
and indeed our numerical analysis agrees with the perturbative analysis even for N'^/q ~ A^. 

The Hamiltonian is made up of the kinetic and potential energy. The potential energy 
actually contributes at order l/g^, which we will neglect. The kinetic energy contains terms of 
order q. These are in principle the dominant terms for large q. They have a constant energy 
density, that is the energy associated to each state| n){n \ . In total the trace gives 

HkinXq) = 2~Nplq. (4.40) 

At order zero, the kinetic energy density is linear (in n), but the total after the trace is zero: 

For the order 1/q contribution we need w{nf' to order l/q^ . 

/ / ^^2 2 Vo/ , 2pi \ 1 / A V ^ 2piV ^ ( AGi (n) 



q \ X J q^ \A \ X J " V 2 

The kinetic energy then is 



(4.42) 



HMnxm = (^Yl - ^ '^'''^') ' ') ^^-^^^ 

which gives 

2- fX^N \ 

Hkin.{l/q) = y (^^(^' - 1) + 2poP2A^j . (4.44) 

Hence the energy up to the 1/g is 

H 2^ 2ENplq + -=- [-^{N^ " 1) + 2poP2N^ . (4.45) 
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The two solutions (for po) behave quite differently for po = X the energy diverges as q becomes 
large while for po = it vanishes. 

Hp,=x ^ 2~NX^q + ^ ( ^(^' - 1) + 2A/?2A^j (4.46) 

^ -EX^NjN^ - 1) 

^''°=° 24^^ ^^-^^^ 

These solutions correspond to an annulus of large radius and small (relatively) thickness that 
oscillate in time. We can see for the operator of the radius squared (which is diagonal), from 
equation (^^), the dominant part of the coefficient (hence eigenvalue of R^) is proportional to 
2q for every state, in the large q limit. The next dominant term is proportional to the square 
root of q and is an oscillatory term (the frequency for po = X \s of order q^ and it is weakly 
n dependent, and for po = it is of order q^^ ). The next important term is proportional 
to 2n + 1. As g S> y/q ^ 1 the radius is very large and oscillates with an amplitude that is 
relatively much smaller than the radius. Thus from afar the droplet looks like an thin annulus 
that undergoes a nontrivial oscillation. 

4.1.2 Small quasi-hole solution 

The small quasi-hole solution, for (g <C 1), can be obtained by a perturbation on the solution 
found in |12|. Take D at first order perturbation to have the form 

D = D' + Dp, Do = D'q + Dop, ^ = ^' + (4.48) 

where the primed variables are solutions found in [|l^] and the variables subscripted p are the 
perturbations. To first order the Gauss law gives 

= [Dp, [D'„ D't]] + [D', [Dop, D'^] + [D', [D'„ Z)t]] + [Z)t , ^d',, D']] 

+ [D'\ [Dop,D']] + [D'\ [D'o,Dp]] + X{[Dp,D'^] + [D',dI]) + + ^'^-t (4.49) 

while the Ampere law gives 

[Dop, [D'o,D']] + [D'o, [Dop,D']] + [D'^, [D'o,Dp]] + [Dp, [D',D'^]] 
+ [D', [Dp, D'^]] + [D', [D', Z)t]] = X[Dop, D'] + X[D'o, Dp] (4.50) 

and the constraint on ^' gives 

Dop^' + D'o^p = 0. (4.51) 
We assume the perturbation takes the form 

Dp = Vge"^ I N - 1)(0 I, Dop = 0, ^p = 0. (4.52) 



The ansatz in |12] has the form (gauge equivalent to) 



N-2 



Z?' = ^ 7GHe™(")* I n)(n + l I, Ao = (4.53) 



n=0 



-12- 



where G{n) and w{n) are explicitly calculated (numerically) in |12|. Then the Gauss law is 
automatically satisfied using 



N-2 



N-1){0 I, f(n) I n+l){n 



n=0 

The Ampere law becomes 



N-2 



0)(7V-1|, J^/(n)|n)(n+l 



?i=0 



0. 



(4.54) 



(4.55) 



with [D', D'"^] given by [4.9| . This implies the same condition as before (removing an overall 
factor Dp) ( [4.16 ) however only for n = A'^ — 1 



p^Dp + (G(0) + G{N - 2))Dp = XpDp. 
For other values of n the condition is automatically satisfied. Then p is given by 



P 



A ± - 4(G(0) + G{N - 2)) 



(4.56) 



(4.57) 



As this expression involves a square root, we see that there is no solution for small A, as we 
will confirm below, in the numerical analysis. This is a gap in the spectrum of allowed q for 
small A. The change in energy is given by 

= -2'ETt{[D'^,dI\[D'^,D'] + [D'^,D'^\[D'^,Dp] + [D'^, dI\[D'^, Dp]) (4.58) 

which gives 

= -2ETr{-pDl[D'o, D'] + p[D'q, D'^]Dp - p^DlDp) (4.59) 
and using that the two first terms vanish, yields 

AH = 2Ep^q. (4.60) 

4.2 Numerical analysis of the quasi-hole solution 

In this section, we will find numerical solutions of equation ( 4.16| ). As our solutions are 
built upon the solutions found in |12|, we first give the numerical analysis of the equations 
considered there. We take this opportunity to correct certain errors that have appeared in 
Equation (5.6) in is incorrect, the kinetic energy is not symmetric about M = N/2. 
The correct equation is: 



/ N-2 \ 

T = 2H { ^ n„ - UnV^Un] +g\N^-N- 2NM) 

\ n=0 / 



(4.61) 



Correspondingly, figure (5) in |12| is also not correct. In addition, there seems to be an incon- 
sistency between figure (5) and figure (7) in |12|, the kinetic energy appears to be greater than 
the total energy. We give corrected figures, figure || and figure ^ here (not for the same values 
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9=0.0625 
9=0.075 
9=0.0S75 



0^ 



^^^*^^' 



9=0.05 
9=0.0625 
9=0.075 
9=0.0S75 



Figure 1: Kinetic energy as a function of M 
for N=100 for (corrected) solution in [|l2| 



Figure 2: Potential energy as a function of 
M for N=100 for solution in [Fl 



of the parameters). These are representative solutions for the system of equations studied 
in Our further numerical analysis concerns solutions built upon these. We ultimately 



use the Newton method |14] to solve the difference equations. Finite difference equations on 
a finite set of variables/parameters and boundary conditions, after iterating recursively and 
removing dependent variables, simply become a system of extremely complicated algebraic 
equations in a much reduced set of variables. In our case, we reduce the system to two equa- 
tions in three variables, which we take to be q, p and G{N — 2). We solve this resultant system 
by the Newton method. We can constrain our search a little by noting that q + G{n) > 
since it appears under a square root in all the expressions for D. For the case n = A*" — 1, if 



we add 2g to the equation ( ^.16 ) then we get the constraint 

"a- VA2 + 8g A + yFT 



(4.62) 



Numerically, we easily find the two solutions that we have determined analytically, for 
q ^ 1. We see in figure (||) that the kinetic energy either diverges as q or vanishes as 1/g 
as g ^ oo, depending on which branch we consider, exactly as we have seen in the section 
(fl.l.ll) . For the potential energy, we give a graph as a function of the value of G{N — 2). The 
preceding two branches of solutions, are found in figure (§), on the bottom part of the curve, 
symmetrically on either side of the point G{N — 2) = 0. This dependence on G{N — 2) can 
be inferred from the analytic solution in section ( [4.1 .11) . We can see in figure (^) that there 
exists a region of transition between the two branches of solutions. This transition region 
seems to imply a lower limit on the permitted value of q for given A. The corresponding 
transition region in figure (|^ is the upper part of the solution between the two peaks. Further 
numerical analysis varying A seems to indicate that the transition region actually extends all 
the way to g = for sufficiently large A. The q = limit corresponds to the case studied 



in ||12|, where it was found that there are N solutions for the noncommutative droplet. We 
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-0.25 -0.2 -0.15 -0.1 -0.05 0.05 0.1 0.15 0.2 0.25 

G(N-2) 



Figure 3: Kinetic energy as a function of q 
for N = 100 and for 4 values oi g = A/2 in 
units of S. 



Figure 4: Potential energy as a function of 
G{N - 2) for N = 100 and g = A/2 in unit of 



find that we recover this multiplicity in the numerical analysis. Specifically the curve of the 
potential energy, figure (^) and figure (||) obtains peaks in the transition region, each of 
which extends down to the solution found in [l^ for different choices of ^, as seen in figures 
(0) and (P). (One has to be careful in this limit for oo \N — 1) . One recovers the solution of 
for {G{N — 2 — M) + q) and not simply q — > 0.) The total energy is always dominated 
by the kinetic term and the range of the total energy is from to oo, however the interesting 
structure in the potential energy implies that its contribution to the specific heat could be 
most important. Comparing our solution to the one found in |12], figures (|l|) and (^), we find 





Figure 5: Potential energy as a function of Figure 6: Potential energy as a function of 
G(A^ - 2) for g = 1 G(iV-2) for5= 1.5 

that the potential energy, figure @, of our solution is at least one order of magnitude smaller 



than that found in ||12| for small A and will approach the potential energy of the solutions in 
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Figure 7: Kinetic energy as a function of q Figure 8: Kinetic energy as a function of q 
for iV = 3 for g — 1.5 



for large A (and small q). The kinetic energy is, however, quite different. In our solution, 
the kinetic energy can take any value since g is a free parameter, while the solution in |12] 
it is constrained to a discrete set of values. For the values of A used in our figures, the 



kinetic energy of the solution in |12| is comparable to the kinetic energy of our solution in the 
transition region. 
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Figure 9: Kinetic density energy for = 100 Figure 10: Potential density energy for N = 
and g = 0.0875 100 and g = 0.0875 

The detail of the transition region is found in the figures (^), ( [To| ) and (^) which describes 
the right half of the transition region {G{N — 2) > 0) in the figure (|^. The potential energy 
density is seen to translate as we vary G{N — 2). This translation continues for G{N — 2) < 
symmetrically with respect to the behavior for G{N — 2) > 0. The potential energy density 
is found to concentrate in the bulk of the droplet, away from the boundaries in the middle of 
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the transition region. During the transition the kinetic density evolves from a hnear density 
{po = A) as in equation ( [4.41 ) to a quadratic density {po = 0) as in equation ( [4.43|) . 

Away from the transition region (i.e. for g ^> 1, we do not give a figure since this region can 
be computed analytically), the form of the potential energy density and the G(n)'s does not 
change markedly, apart from their overall size. The potential energy density is concentrated 
around the boundary. 
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Figure 11: G(n) for TV = 100 and g = 0.0875 Figure 12: Kinetic energy as a function of q 

for g = 0.05 

Finally, in the figure (p^, we show the solutions for various values of N. We see that 
the minimum value of q seems to be growing linearly with A''. Thus, assuming that the trend 
continues, as — > oo then qmin. — > oo, thus the large g-hole solutions would not occur in the 
infinite plane, at least for small A. 

5. Ground state 



The solutions that we have found do not easily allow us to identify the ground state. For 
the infinite case, the ground state is given by the solution in terms of simple annihilation and 
creation operators, |§, 11|, 

a, a^, [a, = 1. (5.1) 

This solution corresponds to a static, quiescent fiuid. However the value of [D, is a non- 
zero constant, which should be considered as the zero point of energy. The solutions in the 
finite droplet approach this state arbitrarily closely in the limit A^ — > oo and for large A. We 
show below in figures ( [T3| ) and ([l^, the plot of G{n) forM = N — l and with A = 2 and the 
corresponding vorticity of the fluid. If we subtract out the constant background due to the 
noncommutativity it is evident that the fluid is almost everywhere quiescent with only net 
vorticity imposed near the boundary. We take these states as the ground states. We also show 
the corresponding kinetic and potential energy in figures ([l5| ) and (16). The kinetic energy T 
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Figure 13: G(n) for M = TV - 1 with N = Figure 14: Vorticity ("magnetic" field B ) for 
100 M = N -1 with TV = 100 (figure cut-off at 

B > -10) 



is concentrated at the boundary, as is the potential energy V if we subtract off the constant 
zero point energy. 
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Figure 15: Kinetic energy density for M 
N - 1 with N = 100 



Figure 16: Potential energy density for M — 
N - 1 with N = 100 (figure cut-off at ^ < 5) 



The quasi-hole solutions that we have found are deformations of these ground states. 
Comparing to the case of a quasi-hole with small we see from equations ( [1.59| ) and ( [4.60| ) 
that the energy density perturbation is solely in the kinetic energy. It is also localized at the 



origin as is evident from the form of the perturbation equation (16). The change in the energy 
is linear in q. The solutions with large q are not perturbative and should be rightly considered 
as solitons. The question of the stability of our identified ground states against creation of 
solitons with large q (the branch with decreasing enenrgy, see figure (^)) is beyond the scope 
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of this paper. 

The rotational excitations that we have found presumably give the Landau levels, upon 
quantization. These are extremely widely separated in energy for large external magnetic field 
and hence will decouple from small excitations. The ground state will correspond to the lowest 
Landau level. The degeneracy of a given Landau level is obtained from the total number of 
particles available and the physical space in which these particles are constrained. The Landau 
levels will presumably play an important role in describing the transitions between different 
Hall plateaux. 



6. Conclusion 



In this article we have found a rich spectrum of the excitations of the noncommutative droplet 
with an action based on the Chern-simons term and an additional Maxwell term. The existence 
of the Maxwell term seems to have profound implications on the spectrum of fundamental 
oscillations. We have found two kind of solutions, the rotational excitation and the quasi-hole 
solution. 

The rotational excitations will exist in the finite or infinite case above any solution of 
the equations of motion. The rotational frequency is exactly the cyclotron frequency for the 
electrons in the model. The rotational excitation correspond to the energy of electrons 
moving in the magnetic field when the amplitude of the rotational excitation is fixed to the 
radius of the cyclotron motion. With the quantization of these oscillations we should recover 
the familiar Landau levels. 

The quasi-holes solutions seem to exist only for the finite droplet. The large quasi-holes 
correspond to a thin annulus that undergoes nontrivial oscillatory vibration. They give a 
continuous band of solutions built above each soliton solution of the field equations. They 
also show a gap in their spectrum as a function of the "charge" q for small enough A. 

We expect that this spectrum of excitations will give rise to a complex phenomenolgy 
which will allow us to describe transitions in two dimensional magnetohydrodynamics, even 
the quantum Hall system. 
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